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ABSTRACT 


Amplitudes  of  first  higher  mode  Love  waves  from  seismic 
events  in  North  America  were  measured  on  the  North  American  continent  and 
corrected  for  attenuation,  structure,  and  distance.  The  corrected  amplitudes 
were  fit  to  a source  model  to  estimate  the  source  depth.  The  high  value  of  the 
attenuation  coefficient  and  lack  of  a good  seismic  model  prevented  good  agree- 
ment between  known  and  calculated  depths. 

In  order  to  estimate  the  effect  of  the  seismic  source  structure, 
the  ratio  of  first  higher  to  fundamental  Love  wave  amplitudes  was  found  ex- 
perimentally and  compared  to  the  theoretical  prediction.  Agreement  was 
better  than  for  the  spectral  fitting  mediod  but  was  still  short  of  being  satis- 
factory. 
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clusions presented  are  those  of  the  authors  and  should  not  be  interpreted  as 
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the  Advanced  Research  Projects  Agency,  the  Air  Force  Technical  Applications 
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INTRODUCTION 


It  is  well  known  that  earthquakes  and  underground  nuclear  ex- 
plosions radiate  energy  in  the  form  of  surface  waves  which  propagate  at  velo- 
cities which  are  the  roots  of  equations  relating  frequency,  seismic  velocity, 
and  propagation  velocity  (Harkrider,  1970).  There  are  an  infinite  number  of 
roots  to  this  equation  both  for  Love  and  Rayleigh  waves,  and  in  principle 
every  seismic  event  is  accompanied  by  energy  propagating  in  each  of  the 
modes  corresponding  to  these  roots.  In  practice,  only  the  fundamental  mode 
(i.  e.  , Rq,  Lq)  with  lowest  velocity  is  excited  with  large  enough  amplitude  to 
be  useful  for  most  detection  and  discrimination  problems.  This  report  de- 
scribes work  attempting  to  obtain  information  about  the  source  from  hitherto 
unused  higher  mode  surface  waves. 

Two  uses  for  higher  modes  are  examined  here,  both  employing 
the  first  higher  mode  Love  wave  (denoted  This  mode  was  chosen  because 

its  amplitude  is  expected  to  be  larger  than  that  of  any  other  higher  mode  Love 
wave,  and  because  it  varies  more  in  amplitude  with  period  than  does  the  first 
higher  mode  Rayleigh  wave,  the  other  possible  candidate  for  this  study. 

A number  of  studies  (Turnbull,  et  al. , 1973;  Turnbull,  1976) 
have  described  a procedure  for  fitting  measured  seismic  spectral  amplitudes 
to  theoretical  source  spectra.  Seismograms  recorded  at  stations  with  good 
azimuthal  distribution  about  the  source  (if  possible)  are  narrowband  filtered 
to  estimate  the  Fourier  component  attributable  at  that  period  to  the  arrival  of 
interest  (Alexander,  1963).  These  spectral  amplitudes  are  corrected  for  the 
earth's  attenuation,  structure,  and  for  distance,  and  the  mean  square  differ- 
ence between  them  and  those  predicted  by  a model  is  minimized  as  a function 
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of  the  model  parameters.  The  parameter  most  sensitive  to  this  fitting  pro- 
cess is  the  source  depth,  so  that  this  method  can  be  used  to  discriminate  be- 
tween deep  and  shallow  events. 

The  radiation  pattern  for  higher  mode  surface  waves  is  simi- 
lar to  that  of  the  fundamental  mode,  since  the  various  modes  are  characteris- 
tic of  the  medium  rather  than  the  source  (Turnbull,  1976).  The  amplitude 
response  of  the  earth  is  different  for  different  modes,  due  to  their  different 
wavenumbers,  and  this  response  can  be  calculated  for  a given  earth  model. 
Finally,  the  attenuation  of  the  earth  may  be  different  for  various  modes,  and 
must  be  found  by  experiment.  Then  higher  mode  amplitudes  may  be  used  to 
calculate  source  depth  in  the  same  way  as  are  fundamental  mode  amplitudes. 

It  was  attractive  to  look  for  another  discriminant  based  on 
higher  modes,  and  the  amplitude  ratio  was  investigated  for  this  pur- 

pose. This  ratio  was  calculated  for  source  models  at  various  orientations 
for  a range  of  periods  from  6 to  15  seconds,  for  several  earth  models.  Under 
the  appropriate  models  the  ratio's  level  changed  with  orientation  but  its  gen- 
eral features  as  a function  of  period  did  not,  and  the  ratio  displayed  features 
below  15  seconds  which  occur  at  lower  periods  for  lower  velocity  surface 
layers.  Measurement  of  this  ratio  should  then  give  an  indication  of  the  seis- 
mic velocity  at  the  source. 
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SECTION  II 

EXPERIMENTAL  t=»ROCEDURE 


The  spectral  fitting  method  outlined  in  the  previous  section 
was  applied  to  data  collected  by  instruments  and  events  from  the  same  con- 
tinent, in  order  to  simulate  the  counter  evasion  problem  as  closely  as  possible. 
The  parameter  determined  by  this  method  is  the  source  depth,  which  is  an 
important  discriminant  between  underground  nuclear  explosions  and  earth- 
quakes. This  experiment  is  described  in  Subsection  A. 

A second  use  for  the  same  kind  of  data  is  presented  in  Subsec- 
tion B on  the  ratio  L^/  L^.  This  ratio  is  expected  to  yield  information  about 
the  source  structure  given  that  the  event  is  shallow.  This  information  is 
supplementary  to  that  required  for  discrimination. 

A.  SPECTRAL  FITTING 

Data  for  this  experiment  were  collected  by  the  Special  Data 
Collection  System  stations.  The  stations  are  listed  in  Table  11-2,  which 
also  gives  their  locations.  Four  events  in  North  America  were  examined 
and  are  listed  in  Table  II- 1.  Another  event  on  the  Burma-India  border  and 
given  in  Table  II- 1 was  also  examined  using  data  recorded  by  the  Seismic 
Research  Observatory  at  Mashhad,  Iran;  the  Alaskan  Long  Period  Array; 
and  at  NORSAR. 

First,  transverse  surface  waves  from  the  Nevada  Test  Site 
(NTS)  explosion  TYBO  were  examined  for  higher  mode  Love  waves,  using 
both  narrowband  filters  and  a maximum  entropy  filter.  It  was  hoped  that 
the  increased  resolution  of  the  maximum  entropy  filter  as  compared  to 
standard  Fourier  techniques  would  lead  to  better  estimates  of  amplitudes. 
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TABLE  II- 1 

EARTHQUAKE  PARAMETERS 


Event 

Date 

Coordinates 

M 

s 

TYBO 

14:00:02 

37.3  N 

4 4 

14  May  1975 

116. 5 W 

MAST 

13:00:01 

37.  3 N 

5 5 

19  June  1975 

116. 3 W 

California 

01:38:51.4 

34.  6 N 

Earthquake 

1 June  1975 

116.4  W 

4.  9 

Wyoming 

18:54:16.5 

44.  9 N 

e 

Earthquake 

30  June  1975 

110.  3 W 

mm 

Burma-India 

05:08:02 

25  N 

Border 

Earthquake 

12  Dec  1975 

96  W 

Hi 

r. 


TABLE  n-2 

STATION  DESCRIPTIONS 


Site 

Location 

Coordinates 

CPSO 

McMinnville, 

38.6  N 

Tennessee 

34.6  W 

FN-WV 

Franklin, 

38.5  N 

West  Virginia 

79.5  W 

HN-ME 

Houlton, 

46.2  N 

Maine 

68.0  W 

WH2YK 

White  Horse, 
Yukon 

60.7  N 
135.0  W 
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However,  it  was  fouad  very  difficult  to  identify  higher  modes  given  only  the 
power  as  a function  of  time  and  period,  which  is  the  form  of  the  maximum 
entropy  filter's  output.  Therefore,  this  approach  was  abandoned  in  favor  of 
narrowband  filtering. 

Next  the  fundaurnental  mode  energy  was  identified  on  traces  fil- 
tered at  periods  between  6 and  35  seconds.  The  apparent  velocity  of  each 
packet  of  energy  so  chosen  was  found  by  dividing  the  epicentral  distance  by 
the  travel  time,  and  plotted  as  a function  of  period  in  Figure  Il-l. 

Also  shown  in  this  figure  are  the  dispersion  curves  for  the 
Hamilton -Healy  and  Gutenberg -Bullen  earth  models.  While  neither  model 
fits  the  data  well,  the  shape  of  the  Hamilton-Healy  curve  is  closer  than  that 
of  the  Gutenberg-Bullen  curve.  For  this  reason,  and  because  it  was  developed 
for  NTS,  the  Hamilton-Healy  model  was  adopted  for  this  study. 

Filtered  seismograms  were  then  examined  for  pulses  of  energy 
arriving  near  the  times  predicted  by  the  Hamilton-Healy  model  for  higher 
modes.  Such  pulses,  if  they  were  clear  and  distinct,  were  identified  as 
mode  surface  waves.  Often,  due  to  the  low  amplitude  of  the  1.^  mode,  a 
clear  identification  was  not  possible.  In  other  cases  the  and  the  pre- 
sumed Lj  energy  overlapped,  making  the  amplitude  determination  dif- 
ficult. In  these  cases  the  identification  and  measurement  of  Lj  amplitudes 
became  to  some  degree  subjective,  and  this  contributed  to  the  experimental 
difficulties  in  this  study. 

Examples  of  filtered  seismograms  are  shown  in  Figure  11-2. 
The  upper  trace  is  the  explosion  TYBO  as  recorded  at  CPSO.  The  outputs  of 
filters  centered  at  10  and  6 seconds  are  shown  in  the  next  two  traces.  Group 
velocities  of  presumed  arrivals  are  indicated  on  the  figure.  While  the  iden- 
tification of  Lq  and  is  fairly  unambiguous  at  6 seconds,  the  presence 
of  Lj  can  only  be  inferred  at  10  seconds  by  the  asymmetrical  shape  of  the 
Lg  energy#  *nd  its  amplitude  is  highly  uncertain. 
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DISPERSION  RELATIONS 


Raw  Trace 


Those  Lj  group  velocities  which  could  be  determiaed  with 
some  certaiaty  are  plotted  versus  period  in  Figure  11>3,  along  with  the  dis- 
persion curves  for  L>j  under  the  Gutenberg -Bullen  model  and  the  Hamilton- 
Healy  model.  Neither  model  fits  the  data  better  than  the  other  over  the  wide 
range  of  periods,  but  the  shape  of  the  Hamilton-Healy  curve  over  the  period 
range  of  interest,  from  6 to  5 seconds,  is  closer  to  the  data  than  that  of  the 
Gutenberg -Bullen  curve.  This  provides  some  additional  support  for  the 
choice  of  the  Hamilton-Healy  model  to  represent  die  propagation  paths  ex- 
amined here. 

Stations  CPSO,  FN-¥TV,  and  HN-ME  are  roughly  in  line  with 
NTS,  as  required  for  the  method  of  determining  the  energy  attenuation  coef- 
ficient of  seismic  waves  described  by  Sun  (1976),  which  assumes  that  the 

-keX 

energy  in  a group  of  waves  varies  with  distance  x as  e Then  a plot 

against  x of  the  quantity 

2 In  A -f  In  sin  ^ 

where  the  term  involving  the  epicentral  distance  A in  degrees  corrects  for 

the  curvature  of  the  earth,  has  slope  -k  Such  a plot  for  L,  amplitudes 

from  TYBO  for  periods  of  6,  7,  8,  9,  and  10  seconds,  where  available,  is 

shown  in  Figure  11-4.  There  is  no  reason  to  believe  from  the  data  that  the 

slope  is  a function  of  period,  so  the  best  straight  line  through  all  the  points  ^ 

above  7 seconds  was  found.  The  slope  of  this  line  is  an  average  value  of  -k^  > 

although  its  intercept  is  difficult  to  interpret. 

The  best  least  mean  squares  fit  to  these  data  gave  k^=  4.  8 x 
10  ^/km.  It  is  difficult  to  assign  a probable  error  to  this  number,  because 
methods  for  computing  standard  errors  assume  that  fluctuations  in  the  data 
are  random,  and  the  dati;  of  Figure  11-4  appear  to  contain  systematic  errors. 

As  a rough  estimate  we  may  suppose  that  the  true  slope  lies  somewhere  be- 
tween the  value  found  for  the  path  CPSO  - HN-ME,  where  k^s  4.  7 x 10~^, 
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FIGURE  II-3 

dispersion  relations 


ENERGY  ATTENUATION 


aad  the  path  FN-MTV  - HN-ME,  where  k * 6.  0 x 10  . Then  our  eetimate 

of  is  4.  8 y X lO'^/km.  This  value  is  au  order  of  magnitude  higher 

than  those  found  by  other  workers  (Tryggvason,  1965),  but  is  consistent  with 
the  fact  that  the  group  velocity  for  this  model  is  relatively  low.  The  Q , or 
quality  factor,  corresponding  to  this  k is  about  20. 

This  large  value  of  k^  is  responsible  for  a great  deal  of  dif- 
ficulty in  using  amplitudes  in  North  America.  It  enters  the  expression 
for  the  amplitude  corrected  for  attenuation  through: 


Y s A e 


k^x/2 


(U-1) 


where  Y is  the  corrected  amplitude,  A is  the  measured  amplitude,  and  x 


is  the  distance.  The  fractional  error  in  Y due  to  errors  in  k is 

e 

k X dk 

dY  e e 


(n-2) 


and  for  the  value  of  k found  here  is  about  5 for  typical  epicentral  dis- 

tances.  Consequently,  fractional  errors  in  the  corrected  amplitude  will  be 
about  five  times  those  in  the  attenuation  coefficient.  Errors  in  the  attenua- 
tion coefficient,  in  turn,  are  expected  to  be  large  due  to  the  uncertainty  in 
Lj  amplitudes. 

Furthermore,  the  value  quoted  for  k^  is  valid  only  for  the 
path  where  it  was  measured,  and  can  only  be  assumed  to  be  valid  for  paths 
to  other  stations  and  from  the  source  to  the  closest  station.  The  use  of  such 
other  stations  is  essential  to  take  advantage  of  the  information  contained  in 
the  radiation  pattern  of  the  seismic  energy. 

This  problem  is  not  a serious  one  for  earth  structures  where 
k^  is  small,  since  the  correction  factor  given  in  Equation  (II-l)  is  nearly 
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Accepting  the  value  of  k quoted,  amplitudes  from  CPSO  and 
WH2YK  were  corrected  for  instrument  response,  attenuation,  and  distance 
over  the  range  of  6 to  10  seconds,  where  energy  from  TYBO  was  observ> 
able.  The  corrected  spectra  were  fit  to  a double  couple  source  model  at  var- 
ious strikes,  slips,  dips,  and  at  various  depths  in  a Hamilton- He aly  earth 
structure,  as  described  by  Turnbull  (1976).  The  model  with  the  smallest 
mean  square  deviation  from  the  observations  at  each  depth  was  chosen,  and 
this  mean  square  deviation  plotted  as  a function  of  depth  in  Figure  11-5. 

Good  agreement  between  model  and  observations  should  be  de- 
tectable by  a minimum  in  Figure  11-5,  but  no  such  dip  is  visible  there.  The 
model  is  nearly  an  equally  bad  (or  good)  fit  to  the  data  at  every  depth.  Sim- 
ilar results  were  obtained  when  amplitudes  from  the  NTS  event  MAST 
and  the  California  earthquake  of  1 July  1975,  were  fit  to  the  model. 

This  lack  of  depth  discrimination  in  the  method  is  due  to  two 
causes.  First,  the  large  value  of  k found  above  means  that  L,  will  have 
a low  signal-to-noise  ratio  at  teleseismic  distances,  and  this  makes  it  dif- 
ficult to  measure  k accurately.  Consequently,  large  errors  will  appear  in 
the  corrected  amplitudes,  which  will  then  fit  no  model  well. 

Another  possible  source  of  difficulty  is  that  the  model  may  not 
contain  enough  details  of  the  earth  structure  at  the  source  region.  However, 
such  fine  structure  would  presumably  lead  to  a sharpening  of  any  minimum 
in  the  residuals  of  Figure  11-6,  udiereas  no  minimum  was  found.  Consequent- 
ly, we  conclude  that  problems  with  the  data  are  great  enough  that  the  use  of 
such  fine  structure  is  unjustified. 

An  event  on  the  Burma-lndia  border  was  examined  at  Mashhad 
(Iran),  ALP  A,  and  NORSAR  in  the  hope  that  the  paths  to  these  stations  would 
display  less  attenuation  than  paths  across  North  America.  Unfortunately,  no 
Lj  energy  was  observed  in  the  filtered  seismograms.  Presumably,  the  rela- 
tively complicated  travel  paths  obscured  higher  modes  at  these  stations. 
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Lj/Lp  RATIO 


Use  of  the  ratio  L^/L^  should  eliminate  problems  associated 

with  inaccurate  determination  of  k.  if  k.  is  the  same  for  and  Li^.  The 

e e 10 

data  of  Figure  U-6  for  TYBO  JL  amplitudes  lead  to  a value  of  k above  7 

-3  ® * 

seconds  of  4.  9 x 10  /km.  The  data  are  somewhat  better  clustered  about 

the  regression  line  than  those  for  L. , and  we  estimate  the  range  of  k to  be 

-3  -3  * ® 

from  4.  8 X 10  to  5.7  x 10  per  kilometer.  Consequently  we  conclude  that 

the  value  of  k^  is  the  same  for  fundamental  and  first  higher  mode  Love  waves 

for  this  transmission  path. 

The  ratio  L^/L^  was  calculated  according  to  a double  couple 
model  at  0. 5 km  depth  for  a number  of  fault  plane  orientations  for  both  the 
Gutenberg>Bullen  and  Hamilton-Healy  models.  The  absolute  level  of  the  ratio 
between  6 and  15  seconds  depended  on  fault  plane  orientation,  but  the  location 
of  a peak  at  7 seconds  for  the  Hamilton-Healy  model  and  a shoulder  near  9 
seconds  for  the  Gutenberg-Bullen  model  was  independent  of  orientation.  This 
ratio  for  two  fault  plane  orientations  is  shown  in  Figure  U-7. 

Experimental  tests  of  this  method  were  undertaken  by  examin- 
ing this  ratio  for  the  events  TYBO  and  the  Wyoming  earthquake  measured  at 
stations  CPSO  and  WH2YK.  The  results  are  plotted  in  Figures  II-8  and  II- 9. 
respectively.  The  results  for  MAST  are  similar  to  those  for  TYBO. 

For  both  explosions  there  is  a peak  in  L^/L^  at  7 seconds 
period  as  predicted  by  the  theory.  For  the  Wyoming  earthquake  the  results 
are  less  clear.  The  absence  of  energy  below  9 seconds,  and  its  pre- 
sence above  15  seconds,  are  in  contradiction  to  the  model  prediction.  Con- 
sequently, while  these  results  are  more  encouraging  than  those  of  the  spec- 
tral estimation  method,  they  are  not  conclusive  proof  that  the  shapes  of  the 
curves  in  Figures  11-8  and  11-9  are  controlled  by  the  structure  at  the  source. 
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ENERGY  ATTENUATION 
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FIGURE  U-7 

Lj^/L  FOR  TWO  EARTH  MODELS 
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FIGURE  n-8 
Lj/Lq  for  TYBO 
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Good  Bigaal-to-aoise  ratio  and  a travel  path  within  North 
America  from  a source  region  with  high  seismic  velocity  are  required  for 
a fair  test  of  the  method. 
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CONCLUSIONS 


Surface  waves  from  a number  of  events  have  been  examined 
for  the  presence  of  higher  mode  energy.  Amplitudes  of  the  first  higher  mode 
Love  wave  have  been  measured,  corrected  for  distance,  instrument  response, 
and  attenuation,  and  fit  to  a model  of  the  fault  plane  and  propagation  path.  The 
results  suggested  that  this  mode,  the  most  promising  for  the  purpose,  has  no 
capability  for  depth  discrimination,  at  least  in  North  America.  The  reasons 
for  this  failure  are  the  high  value  of  the  energy  attenuation  coefficient  and  the 
lack  of  a good  source  modeL 

In  order  to  discriminate  between  different  source  structures, 
the  ratio  L^  / L^  was  calculated  as  a function  of  period  both  experimentally 
from  the  observed  data  and  theoretically  from  two  source  models.  Agree- 
ment between  observation  and  theory  was  more  encouraging  than  for  the 
spectral  fitting  method,  but  was  still  short  of  satisfactory. 

It  is  recommended  that  any  further  study  of  the  utility  of  higher 
modes  for  discrimination  be  directed  toward  data  from  a region  of  low  attenua- 
tion coefficient  and  well  known  simple  structure. 
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